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Time-dependent density functional theory (TDDFT) 
is a theoretical framework which, similarly to the ground 
state DFT relies on the one-to-one mapping of the 
density of particles to the external potential [2[. The 
unique density-potential correspondence implies a possi- 
bility to calculate the exact time-dependent density by 
solving Hartrec-likc equations for fictitious noninteract- 
ing Kohn-Sham (KS) particles. This tremendous simpli- 
fication of the problem makes TDDFT one of the most 
popular ab initio approaches for describingquantum dy- 
namics of realistic many-body systems 

Standard TDDFT is formulated for systems of quan- 
tum particles driven by classical electromagnetic fields 
which covers most traditional problems in physics 
and chemistry. However, nowadays the experimental sit- 
uation is rapidly changing. Progress in the fields of cav- 
ity and circuit quantum electrodynamics (QED) opens 
a possibility to study many-electron systems strongly 
interacting with quantum light. Notable examples are 
atoms in optical cavities in cavity- QED 0-0]' or meso- 
scopic systems such as superconducting qubits [iWlol] and 
quantum dots [8 10] in circuit-QED. Recently a strong 
coupling of molecular states to microcavity photons and 
the modification of chemical landscapes by cavity vac- 
uum fields have been reported [Ill4l3j . Obviously, the 
classical treatment of external fields prevents application 
of TDDFT to this new interesting class of problems. 

This paper presents TDDFT for systems of electrons 
strongly coupled to (or driven by) a quantized electro- 
magnetic field (3]. We prove the generalized mapping 
theorems for TDDFT and for time-dependent current 
density functional theory (TDCDFT). In both cases we 
analyze the structure and properties of the exchange cor- 
relation (xc) potentials and discuss possible approxima- 
tion strategics. Finally we make a connection of the 
present theory to TDDFT for open quantum systems. 

Consider a system of N electrons, e. g., an atom or a 
molecule, placed inside a cavity hosting M photon modes. 
In the Schrodinger picture the configuration of the system 



is specified by the positions {xj}^^ of the electrons and 
the set {qa)a=i of photonic coordinates. The full sys- 
tem is described by the wave function >F({xj}, {q a }, t). 
Assuming as usual [15| that the wavelength of relevant 
photon modes is much larger than the size of the elec- 
tronic system we adopt the dipole approximation for the 
electron-photon coupling. The Hamiltonian of the sys- 
tem takes the following form 
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where W Xi - Xj is the electron-electron interaction, u a are 
the frequencies of the photon modes, and X a describes 
coupling to the a-modc. The electron and photon sub- 
systems can be driven externally by the classical vector 
potential A oxt (x,i) and the external "currents" J" xt (t). 
The vector potential describes (in a temporal gauge) 
forces from the ions in atoms and molecules and all other 
possible classical fields. The currents J" Kt (i) allow for an 
external excitation of the cavity modes. The time evolu- 
tion from a given initial state ^({xj}: {la}) is governed 
by the Schrodinger equation 



(2) 



Solution of this equation gives a complete description of 
the system for a fixed configuration of the external fields. 

All DFT-like approaches assume that the state of the 
system is also uniquely determined by a small set of ba- 
sic observables, such the density in TDDFT, the cur- 
rent in TDCDFT, and possibly something else in other 
generalizations of the theory. Below we construct two 
generalizations of TDDFT for electron-photon systems 
described by the Hamiltonian of Eq. fl). 

Let us start with a technically simpler current-based 
theory. Consider the electronic current j(x, t) and expec- 
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tation values Q a (t) of photon coordinates as the basic 
variables. These variables arc defined as follows 



(*l9a|*>, 



(3) 



j = (*|] p (x)|*> - -A ext - V ^<*|g a fi(x)|*), (4) 

fU ' 7TJ. 



where n(x, t) = {ty\h(~x)Y&) is the electronic density, and 

™( x ) = Ej 5 ( x - x j) and Jp( x ) = Ej{ v ^ 5 ( x - x i)l 
are the density and paramagnetic current operators. 

Equations of motion for Q a follow Eqs. © and ©: 



\ a 3(t) + J^ t (t) 



(5) 



where 3(t) = J j(x, t)dx is the space-averaged electronic 
current. Equation © is simply the Maxwell equation for 
the cavity vector potential projected on the a-mode. 

Equations (JT|)— (|4j) determine the wave function ^(t), 
and the basic variables j and Q a as functionals of the ini- 
tial state Vfo, and the external fields A ox t and J^ t . This 
defines a unique map {* , A ext , J" xt } H> j, Q a }. TD- 
CDFT assumes the existence of a unique " inverse map" 
{*o,j,<9a} !->■ {^7 A ext , J^ t }. That is, given the initial 
state and the basic observables one can uniquely recover 
the full wave function and the external fields that gener- 
ate the prescribed dynamics of the basic variables. 

To prove the uniqueness of the inverse map we follow 
the nonlinear Schrddingcr equation (NLSE) approach 16, 
17] . Assume that j(x, t) and Qo.it) are given, and express 
the external fields from Eqs. © and © as follows, 



A ext = ~<* |j P (x) - j|*> - J2 — (*|g Q n(x)|*), (6) 

a 

J c " xt = Qa + U*Q a - X a J. (7) 



This defines the external fields as explicit functionals of 
the observables j(x, t) and Q a (t), and the instantaneous 
state V(t). Substitution of Eqs. © and into the 
Hamiltonian ((T|) turns Eq. ([2]) into the many-body NLSE 
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where H[j, Q a , \&] is an instantaneous functional of ^(t), 
which depends paramctrically on j(x, t) and Q a (t). The 
uniqueness of a solution to Eq. © can be proven eas- 
ily under the usual in TD(C)DFT assumption of t- 
analyticity 

n m m. Assumm s tnat j( x ^) and Qa(t) 

are analytic functions in time, we represent them and 
unknown ^(t), by the Taylor series 



m = E j (fc)t "' = E Q« kk > *(*) = E * 
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After inserting these series into Eq. © one observes that 
all coefficients with k > can be expressed recur- 
sively in terms of j( k \ Qa , and = ^q. The simple 
reason for this is that the right hand side in Eq. © is 



an instantaneous functional of ^(t), while the left hand 
side ~ dt^(t). As the recursion produces the unique 
Taylor series for ^(t), the many-body wave function is a 
unique functional of the initial state and the basic vari- 
ables, VPfv^o Jj Qa]- By substituting this wave function 
into Eq. © we find the functional A cx t[\&o,j, Qa], which 
completes the proof of the TDCDFT mapping theorem. 

The KS system for this theory is constructed as follows. 
Consider a system of N noninteracting particles coupled 
to the photon modes at the mean field level. This system 
is described by a set of N one-particle KS orbitals tpj (x, t) 
which satisfy the following equations 
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where Q a {t) is the solution to Eq. © with 3(t) being 
replaced by the space-average of the KS current density 
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Using the above NLSE argumentation (or the standard 
TDCDFT mapping [l^]) we find that <fj(x, t) are unique 
functionals of the KS current js(x, t) and the KS initial 
state \&o • A comparison of Eqs. © and (fTU]) shows that 
the KS current reproduces the physical current js = j if 
As in Eq. © is defined as As = A oxt + Ah xc , where 



Alixc = „ [E Im (^* V ^) ~ m (*ljp( x )l*) 



5^A a (tf|Ag a Aft(x)|¥) 
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Here A.q a = q a — Qa{t), and An(x) = h(x) — n(x, i) 
are the fluctuation operators for the photonic coordinates 
and the electronic density, respectively. By construction, 
for given initial states and "Fg , the potential Ah xc is 
a gauge invariant universal functional of j and Q a . 

Therefore the current and the photonic coordinates can 
be calculated from a system of Eqs. ©, © describ- 
ing noninteracting fermions driven by a selfconsistent 
field and coupled to a set of classical harmonic oscilla- 
tors. There is a deep reason for explicitly separating the 
"mean-field" part A m f = ^2 a X a Q a of the selfconsistent 
potential in Eq. ©. This part accounts for the net force 
exerted on electrons from the photons. The remaining 
Hxc-part Ajj xc does not produce a global force, 



J [j x (V x A Hxc ) - nd t A Hxc ]dy 



0, 



(12) 



which can be checked directly using Eqs. (|TT|) . ©, and 
©. Apparently both Hartree and xc contributions to 
Ahxc = Ah+A xc satisfy the identity of Eq. (fl~2"|) indepen- 
dently. Equation (|T2"j) is a generalization of the zero-force 
theorem [201 ] for the considered electron-photon system. 
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In this regard A xc is similar to the xc potential in the 
usual TDCDFT for closed purely electronic systems. 

If the electrons are driven by a scalar external poten- 
tial 14 x t(x, t), one can choose the density n(x, t) as the 
basic variable for electronic degrees of freedom. Let us 
identify the photonic basic observables for the electron- 
photon TDDFT. First, we transform of the photon field 
in Eq. (UJ from the velocity to the length gauge flBT ]. 
Then we perform the canonical transformation of pho- 
ton variables, id qa >->• w a p a , q a h- > —iui' a l dp al which 
"exchanges" the photon coordinates and momenta while 
preserving their commutation relations. The system is 
now described by the wave function $({xj}, {p a }, t) that 
is the Fourier transform of ^({x,,}, {q a }, t). The final 
transformed Hamiltonian takes the form 
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where X = J2j=i x i- The structure of Eq. (fl"3| suggests 
that the proper basic variables are the density n(x, t) and 
the expectation values P a (t) of the photon momenta 

n(x,t) = ($|n(x)|$), P a (t) = ($|p a |*). (14) 



Equations of motion for the basic variables read 

P a + u 2 a P a - cj a \ a R = - J c a xt /w Q , (15) 
rnii + VF str + Vf « = v (" vt/ "cxt), (16) 



where R(f) = ($|X|$) = f xn(x, i)<fx is the expecta- 
tion value of the center of mass coordinate, and F str = 
im($| [T + W, j p ] |$) = — VII is the usual electronic stress 
force which is equal to the divergence of the electronic 
stress tensor. The force f Q (x, t) exerted from the a- 
photon mode on electrons is given by the expression 

k(x,f)=A a <$|(w a p„-A a X)n(x)|$). (17) 

Now we are ready to prove the uniqueness of the map 
{&o,n, P a } h-> {$, V ex t, Jext} f fom the initial state and 
the observables to the time-dependent wave function and 
the external fields. The corresponding many-body NLSE 
is constructed using Eqs. (fT5)) and (fT6|) . 

By solving Eqs. (fl"5)) and ([To]) for J£ xt and V ext we get 
the external fields as functionals of the basic variables 
and the instantaneous wave function $(i), 14 x t[$, n] and 
^x t [n, Pa\. Inserting these functionals into Eqs. (fi"3|) 
we obtain a ^-dependent Hamiltonian H[n, P a , <fr] of the 
many-body NLSE that determines the wave function for 
a given initial state $o, the density n(t) and photon mo- 
menta P a (t) . The uniqueness of a solution to this NLSE 
is demonstrated in exactly the same way as for the above 



TDCDFT case, provided the standard i-analyticity con- 
ditions arc fulfilled. This proves the generalized TDDFT 
mapping theorem: the many-body wave function and the 
external fields are the unique functionals of the basic vari- 
ables, n and P a , and the initial state (2lj . 

The KS system can be again constructed explicitly. 
Consider a system of noninteracting particles described 
by N KS orbitals which satisfy the equations 



id t (j>j 



F g + ^(o; a P a -A a R)A a x]^, (18) 



where the second term in the square brackets is the mean- 
field analog of the electron-photon interaction term in 
Eq. ([13"]) . The force balance equation for this system 
takes the form 

mn + VFf tr + V ]T \ a {uJ a P a ~ A Q R)n = V(nW«j), 



where F^ r 



(19) 

im($ s \[f , j p ]\$ s ) = -VLIs is the kinetic 
stress force of noninteracting fcrmions [$ s (t) is the KS 
Slater determinant]. By applying the NLSE arguments 
to Eqs. (fT8"|) - (fT9"|) we conclude that y>j and Vs are unique 
functionals of n(x,i), P a (t), and the KS initial state $q. 
Then, from Eqs. (JTBj! and one finds that the KS 
density reproduces the exact density if Vs is of the form 



V s = Kxt + V, 



Hxc 



ya 



(20) 



where the universal functionals ^^[n, P] and V£.[n,P] 
are defined via the following Sturm-Liouville problems 

V(nV< c ) = V(Ff tr - F str ) = V(vShxc), (21) 
V(nW X c) = VA Q {$|(A Q AX-w Q Ap a )An|$), (22) 

with Ap a = p a — P a (t) and AX = X — R(i) being the 
fluctuation operators for the photon momenta and the 
center of mass coordinate of the electrons. 

Interestingly, in TDDFT the total xc potential is nat- 
urally separated into the usual electronic stress contribu- 
tion and the contributions assigned to each photon mode. 
It is obvious from Eqs. (|2ip and (f2"2"j) that each contribu- 
tion to the total xc potential satisfies the zero-force the- 
orem, 

/ nW Hxc rfx = / riVV™ c dx = 0. The net photon 
force exerted on electrons is fully captured by the mean 
field electron-photon potential in Eq. (fTBf . The zero force 
theorem is a consequence of the harmonic potential theo- 
rem (HPT) [Ulll^li which also holds true here as photons 
form a set of harmonic oscillators coupled bilinearly to 
the electronic center of mass. 

In practice any DFT-type approach requires approxi- 
mations for xc potentials. In the present generalization 
of the theory we succeeded to define the xc potential in 
such a way that it has the same general properties, and 
obeys the same set of constraints as the xc potential in 
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the usual purely electronic TDDFT. This suggests natu- 
ral strategics for constructing approximations. 

(i) The first possibility is a velocity gradient expan- 
sion. At zero level we set Vg. = and take = V^ LDA , 
the xc potential in the standard adiabatic local density 
approximation (ALDA). This seemingly naive approx- 
imation exactly reproduces the correct HPT-type dy- 
namics as for the rigid motion with a uniform velocity 
v = j/n the effect of photons on the density dynamics 
is exhausted by the mean-field contribution. The zero 
level approximation can be viewed as a generalization 
of ALDA. The dynamical corrections should be propor- 
tional to velocity gradients. It should be possible to de- 
rive them perturbativcly in the TDCDFT scheme along 
the lines of the Vignale-Kohn approximation 23J, |24 1 . 

(ii) Probably a more promising strategy is to make a 
connection of the effective KS potential to the many-body 
theory (25[. Beyond the mean field level the electron- 
photon coupling generates a retarded photon-mediated 
interaction between the electrons. The corresponding 
photon propagators will enter the diagrams for the elec- 
tronic self energy as additional interaction lines. The 
new contribution to the self energy can then be con- 
nected to the xc potential via Sham-Schliitter equation 
26l 27 1 . In principle the corresponding xc potential can 
be constructed perturbatively to any desired order in the 
coupling constant f28| . However, already the simplest 
approximation, generated by the exchange-like diagram, 
is expected to capture the important physics. Formally 
this approximation for V" is an analog of the x-only opti- 
mized effective potential 27 , 2!| . Physically it should be 
responsible for the Lamb shift effects and for the sponta- 
neous photon emission in nonequilibrium situations. 

Exploring practical performance of these approxima- 
tions is an interesting direction for the future research. 

If the functionals Vhxc[?"i, Pa] and Vg c [n, P a ] are known, 
the basic variables, n(x, t) and P a (t), can be calculated 
by solving Eqs. (US), (JT5). In general the KS Eq. (JTS} 
should be solved numerically, while Eq. (fTS"]) always ad- 
mits an analytic solution. For example, for the equilib- 

this solution reads 



rium initial state and J^ t 



P a (t) 



sm[u a (t - t')]\ a R(t')dt' . 



(23) 



By substituting Eq. (|23| into Eq. (|18|) we eliminate the 
photon variables and get the KS equation involving only 
the electronic density. Thus we obtain TDDFT for an 
open quantum system - now the KS equation describes in 
a closed form only the electronic part of the full electron- 
photon system. 

Formally Eq. (fl3 ) is a version of the Caldeira-Leggett 
(CL) model |3Cj, |31(. Therefore, as a byproduct we ob- 
tained TDDFT for open systems coupled to the CL bath 
of harmonic oscillators. Let us assume Ohmic spec- 
tral density of the bath, TrJ2 a KKH^ ~ w Q ) = 2t^, 



where i] is the friction constant. In this case the self- 
consistent potential in Eq. (fT8|) reduces to the form 
V e ff = Vhxc + ?7-/VRx. The last, mean field term is exactly 
the potential in the phenomenological dissipativc NLSE 
proposed by Albrecht 32 , 33| . Hence already at zero level 
we recover one of the heuristic theories of quantum dis- 
sipation. Deficiencies of the Albrecht equation should be 
corrected by going beyond the zero level approximation. 

Currently there are several formulations of TDDFT 
for open systems, based on the master equation for the 
density matrix |34l436| , or on the stochastic Schrddinger 
equation 37J, |38[ . At the level of the final KS equations 



our theory is similar to the formulation of Refs. [35 . 
which also allows for the unitary propagation of KS or- 
bitals. The conceptual difference is that in the present 
case both the TDDFT mapping and the approximation 
strategies are universally valid for the cavity situation 
with a few discrete photon modes and for the bath with 
a continuous spectral density. The bath in traced out at 
the very last step after setting up the TDDFT framework 
together with approximations. 

In conclusion, TD(C)DFT for systems strongly cou- 
pled to the cavity photon fields is proposed. We proved 
the corresponding generalizations of the mapping the- 
orems, established the existence of the KS system, and 
suggested a few technically feasible approximation strate- 
gies. In the limit of dense spectrum of photon modes this 
approach naturally leads to TD(C)DFT for open quan- 
tum systems. This work is a step towards ab initio theory 
of various cavity /circuit QED experiments, and practical 
TDDFT for dissipative systems. 
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